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Sensitivity of mixing layers 
to three-dimensional forcing 

By J. J. Riley 1 , P. D. Mourad 1 , R. D. Moser 2 AND M. M. Rogers 2 


1. Introduction 

It is well-known that turbulent mixing layers are dominated by large-scale, fairly 
coherent structures, and that these structures are related to the stability charac- 
teristics of the flow (see, for example, Ho & Huerre 1984, for an excellent review 
of this subject). These facts have led researchers to attempt controlling such flows 
by selectively forcing certain unstable modes, which can in addition have the effect 
of suppressing other modes. For example Oster & Wygnanski (1982) and Ho & 
Huang (1982) found that, by subjecting the mixing layer to oscillatory disturbances 
upstream of the splitter plate at frequencies related to unstable modes, the growth 
of the layer could be significantly inhibited. Similarly it has been found (Husain 
& Hussain 1986) that forcing at certain frequencies and their subharmonics can 
enhance the growth of the mixing layer. 

Much of the work on controlling the mixing layer has relied on forcing two- 
dimensional instabilities. In this study we address the results of forcing three- 
dimensional instabilities. The objectives of our work are twofold: (1) to understand 
how a mixing layer responds to three-dimensional perturbations, and (2) to test 
the validity of an amplitude expansion in predicting the mixing layer development. 
The amplitude expansion could be very useful in understanding and predicting the 
three-dimensional response of the flow to a variety of initial conditions. 

The mixing layer is sensitive to several aspects of the imposed perturbations, 
including their wavelength or frequency (see, e.^., Michalke, 1964), the relative 
phases of the disturbances if more than one are used (e.</., Patnaik, Sherman & 
Corcos 1976), the relative amplitudes of the disturbances (Metcalfe et al 1987), 
and their shape (e.p., eigenfunctions of the Orr-Sommerfeld equation, etc.). The 
approach we would like to take is to systematically explore the sensitivity of the 
mixing layer to various types of three-dimensional perturbations, using both direct 
numerical simulations and amplitude expansions. However, the potential parameter 
space to explore is very large. In this preliminary study we begin by selecting several 
different conditions for simulation and examining the behavior of the computed flow 
fields. In order to test the amplitude expansion technique, we compare previous 
simulation results with predictions from the theory for cases where the simulations 
and theory treat similar initial conditions. In this preliminary study, we limit the 
application of the expansion technique to two-dimensional mixing layers. 
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2. Mathematical description 

Consider a temporally evolving mixing layer with a Cartesian coordinate system 
oriented with x in the mean flow direction, y in the direction of the variation of the 
mean velocity, and z in the transverse direction. 

It is both convenient and natural to consider initial perturbations of the form 

u i (x,y,z) = u i (c t ,y, 1 )e* a * + ^ 

since this is the form of the stability eigenfunctions. The simulations and analysis 
were performed using initial conditions consisting of the sum of a few perturba- 
tions of this form. Consider first the two-dimensional modes, where 7 = 0 . If 
the simulations were initialized at a very low level with broad-banded noise, then 
the first significant disturbance would appear at approximately (aF>0), where ap 
is the wave number of the most unstable mode based upon linear stability theory 
(ap « 0.8892$, where 6 is the initial vorticity thickness, Michalke, 1964). Initializa- 
tion of this mode causes roll up to occur (Patnaik, Sherman & Corcos 1976; Riley 
& Metcalfe, 1980), and the development of higher harmonics (e.^., (2aF,0)). If 
the initial field is broad-banded, then the subharmonic (1/2 of, 0) will also show 
significant growth, both due to its own inherent instability and also to nonlin- 
ear interaction with the fundamental. Vortex pairing will eventually occur. The 
behavior of these two-dimensional disturbances is well-documented by laboratory 
experiments, numerical simulations, and theory. 

Next consider three-dimensional disturbances (7 ^ 0). In laboratory experiments 
used to study the mixing layer downstream of a splitter plate, it is observed that 
the initial primary three-dimensional disturbance consists of streamwise vortices of 
a particular wavelength (approximately 7 f = 1.5of) that appear locked in phase 
(Bernal, 1981; Ho & Huang, 1982). In the current nomenclature, these modes are 
the points (0,7f) in the ( 0 , 7 ) plane. Initial disturbances of this type were con- 
sidered by Pierrehumbert & Widnall (1982), and their growth termed translative 
instability. It was found that, based upon linear stability theory, these translative 
modes are weakly algebraically unstable. However, in the presence of a finite ampli- 
tude, two-dimensional disturbance, the translative modes become highly unstable. 

Another mode that has received some previous attention is the oblique mode at 
(of, 7 f)j which was addressed theoretically by Pierrehumbert & Widnall (1982) 
and Corcos & Lin (1984), and numerically by Metcalfe et al. (1987). The latter 
work showed that initializing these modes along with the fundamental results in 
streamwise disturbances very similar to those observed in the laboratory. This mode 
is unstable according to linear stability theory, but the quoted previous work shows 
that it is much more unstable in the presence of finite amplitude, two-dimensional 
disturbances. 

At this point it is important to mention certain differences between the spatially- 
evolving and temporally-evolving mixing layers. The former is what is usually 
studied in laboratory experiments, and is a more realistic approximation to many 
flows encountered in technology, e.g in mixing regions in reaction chambers. The 



I 


Three’ dimensional forcing of mixing layers 93 

temporally-evolving mixing layer is probably a better model for geophysical situ- 
ations, e.g ., mixing layers in the atmosphere or oceans. It has been studied more 
often numerically because of the significant numerical advantages in avoiding in- 
flow/outflow boundary conditions. Some laboratory experiments have been con- 
ducted for the temporally-evolving case (see, e.y., Thorpe, 1985) using salt-stratified 
water in a tilting tank. 

The temporal mixing layer studied in the tilting tank experiments shows lit- 
tle evidence of the strong streamwise vortices which are commonly observed in 
spatially-evolving mixing layers (Thorpe 1985). This is perhaps due to the fact 
that the (0,7f) mode would not be preferentially excited in these experiments. In 
spatially-evolving mixing layer experiments, streamwise vortices are often present 
at the edge of the splitter plate, and are probably related to features of the particu- 
lar experimental facility (Ho, 1988). Although the temporally-evolving flow exhibits 
the same qualitative two-dimensional features as the spatially-evolving mixing layer, 
e.y., rollup and pairing, it appears that the development of the three-dimensional 
disturbances may be significantly different, due to differences in the nature of the 
disturbances. In the numerical simulations presented in §4, initial disturbances 
which would produce streamwise vortices similar to those observed in spatially de- 
veloping mixing layers were selected. 

To address the problem of the sensitivity of the mixing layer to three-dimensional 
disturbances, we consider the constant density Navier-Stokes equations: 

~ +u- Vu = -Vp+ir^u (2.1) 

V • u = 0 (2.2) 

Here u = (u,v,ti;) is the velocity vector and p is the pressure. The equations have 
been nondimensionalized using the velocity scale A17/2, half the mean velocity 
difference across the layer; the length scale 6 = A U/U y (y = 0) (initial vorticity 
thickness) and the time scale 2S/AU. The Reynolds number R is defined R = 
AUS/2u. 

The initial velocity field consists of an initial mean velocity U plus an initial per- 
turbation u'. U is taken to be an error function, (U(y) = et{(y/Try)) For boundary 
conditions we assume periodicity in the x and z directions, and assume that all 
perturbations decay as |y| —> oo. 

3 Amplitude expansion technique 

We have extended a procedure for theoretical analysis of two-dimensional wave- 
wave interactions in unstable shear flows (Mourad, 1987; Mourad & Brown, 1988) to 
the case of two- and three-dimensional waves, here applied to the temporal mixing 
layer. The details of this derivation and results are in the Appendix. As a first ap- 
plication to the mixing layer, we will consider the classic problem of the interaction 
of a two-dimensional wave with its two-dimensional subharmonic, to determine if 
we can reproduce the phase- dependency of their evolution, as described in Patnaik, 
Sherman & Corcos (1976), Riley & Metcalfe (1980) and Ho and Huerre (1984); see 
also Mansour et al. in this volume. 
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8.1 Two-dimensional amplitude expansions 

We begin by listing the equations governing the temporal evolution for the stream 
function of the two dimensional fundamental, with amplitude and phase A(t) and 
0 A (t), respectively, and of the subharmonic, with amplitude and phase 8 (f) and 
Os(t)> respectively. 

- 7 - 7 - — = oo + ai— + a$A 2 + a 4 B 2 

A at at A 

— — i -^ Ob = 60 + biAe ~ l ( $A ~ 20B ) + 63A 2 + 64-83 

+ b 1 A 2 e- 2i ^ A - 20B) + b s A 2 e-*° A - 2$B) cos(0 A 
+ b 9 A 2 e-* $A - 2$B) sin(0 A - 20 b ) 

As discussed in detail in the Appendix, the rational for these equations is similar 
to that of Lorenz, or of Stuart’s energy method (Ho & Huerre, 1984). Like them, we 
derive Equations (2.3) and (2.4) without benefit of asymptotic analysis, which would 
ordinarily allow for an ordering of the weights of the various physical processes in 
the mixing layer. We forgo asymptotics, because a reasonable small parameter is 
not available for our use in supercritical mixing layers. Instead, the expansions 
are truncated based on the physical processes we wish to include or remove from 
the model. Finally, like Lorenz and Stuart, we hope that the results lend at least 
qualitative insight into the true physical processes. We do, however, include more 
detail in the expansions than these other theories. 

In 2.3 and 2.4, ao and b 0 govern the linear growth of their respective disturbances. 
They are complex constants of the form <r p — u> p t, where cr p is the growth rate and 
u> p is the linear frequency of each of the waves. 

The terms beginning with a\ and 6 j represent the effects of the changes induced by 
each wave on the other. This direct interaction was the focus of the asymptotic study 
of Monkewitz (1988), who used an expansion near the critical Reynolds number 
(% 0) to define a small parameter. However, for supercritical flows, other nonlinear 
effects are also important (see Ho & Huerre 1984). Only a\ and 61 terms are included 
by Monkewitz and others, because they show that when the flow is near critical the 
return-feedback effects discussed below occur on a much longer time scale than the 
initial interaction. 

Next, the term with coefficient a$ represents the affects on the fundamental wave 
of the modified mean flow and the first harmonic (2aj?,0). Both the modified mean 
and the first harmonic are forced by the fundamental. This term generally acts to 
check the unbridled growth of the wave in a supercritical flow. 

The term with a 4 represents a variety of physical effects. These include the effect 
on the fundamental of ( 1 ) the modification to the mean flow induced by the sub- 
harmonic wave; ( 2 ) the distortion of the subharmonic induced by the fundamental 
(as represented by the 6 j term). The energy methods would have a term like a 4 in 
their evolution equations which would not include the second contribution. 


(2.3) 

(2.4) 


— 29b) 
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FIGURE 1. Evolution of the energy in the (A) fundamental (i?( ajPj0 )) and (B) 
subharmonic as predicted by equations (2.3) and (2.4). 

The term beginning with 63 represents the effects on the subharmonic of the 
changes in the mean flow induced by the fundamental, and the 64 term is analo- 
gous to the a3 term. In addition, 64 includes additional effects due to the direct 
interaction of the fundamental and subharmonic, as in the 04 term. 

The more exotic looking terms with bj y b$ and 69 are all due to the impact of the 
direct interaction of (aj?,0) and (ajr/2,0) on the subharmonic. 

S.t Two-dimensional amplitude expansion results 

For this report, we will not consider a detailed discussion of the roles played by 
the various terms in the calculations; though this is of great interest, it is left for 
future work. Here we will verify the effect on the interactions of changing the initial 
phase difference between the fundamental and subharmonic. The importance of 
the relative phase is well documented (Husain & Hussain 1986, Patnaik, Sherman 
& Corcos 1976, and Riley & Metcalfe 1980). Following the definitions used in the 
numerical simulations, the fundamental and subharmonic are in-phase if alternate 
cores of the fundamental modes are aligned with the cores of the subharmonic 
modes. They are out-of-phase if the cores of the subharmonic are between the 
cores of the fundamental. (This is apparently opposite to the definition used in the 
review article of Ho & Huerre, 1984.) An operational definition like that above is 
important, since the linear stability eigenfunctions have an arbitrary phase. 

Figure 1 shows a typical result of integrating the evolution equations (2.3 and 
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FIGURE 2. Contours of the stream function at time 13 as computed using the 
amplitude expansion for disturbances which are in-phase. 


2.4). What is plotted is an approximation to the energy (E^ p ^) in the fundamental 
and subharmonic modes, as given by; 


£(«„0) = A 2 (t) f [«o,(a„0) + »0,(«„0)1 d v 

J — oo 

/ oo 

[ tt o,(a,./2 ) o) + “’o, (or/2,0)] d y 

-oo 


(3.1) 


Note that these are unnormalized energies. The theory is also capable of generating 
a more complete energy integral, with the total nonlinear form U( aF>0 )> rather than 
the structurally linear form u 0> (aFiO) (see the appendix for an explanation of the 
notation). This is a matter for future work. 

The linear growth of each of the waves (the straight lines on this log-linear graph) 
is followed by the weak acceleration of the subharmonic (between nondimension al 
times 13 and 16) due to the effects of the b\ term. The fundamental then equilibrates 
due to the action of a$ y while B keeps growing due to the forcing of A y as well 
as its own ability to extract energy from the mean flow. Once the subharmonic 
amplitude reaches a critical value, the fundamental is reduced, and the subharmonic 
equilibrates due to the action of 64 and 04. (Note that the final state as predicted 
here is suspect, since we have discovered that 64 and a\ are incomplete as currently 
calculated. This will be amended in future work.) This graph is representative of 
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FIGURE 3. Contours of the stream function at time 18 as computed using the 
amplitude expansion for disturbances which are in-phase. 


both the in-phase and out-of-phase calculations, at least qualitatively, and is similar 
to that in Riley & Metcalfe (1980). 

During the linear portion of the evolution, the total stream function field (Figure 
2) shows the characteristic scales associated with the fundamental mode, here at 
time 13 for the in-phase case. The subharmonic is there, with its core at the 
center of the plot underneath that of the fundamental. However, the subharmonic 
makes a negligible contribution to the field at this time. The stream function field 
for the out-of-phase case is the same. Note that this and all the following plots 
have streamwise domain lengths which equal the wavelength of the subharmonic: 
therefore, two fundamental vortices and one subharmonic vortex can fit into the 
domain. 

Once the nonlinear interactions become important (nondimensional time of about 
18), the flow patterns for the in-phase and out-of-phase cases become different. The 
total stream function field for the in-phase case (Figure 3) shows an enhancement 
of the original, central, fundamental vortex, at the expense of the vortices to either 
side, whose torn remnants are at the edge of the graph. This is the ‘tearing’ mode 
discussed in Patnaik, Sherman & Corcos (1976). Note that the mixing region, as 
well as the vortex itself, is smaller here than for the out-of-phase case at the same 
time (Figure 4). There the subharmonic vortex, made up of two paired, fundamental 
vortices, fills the graph, showing an increase in the mixing region. 
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FIGURE 4. Contours of the stream function at time 18 as computed using the 
amplitude expansion for disturbances which are out-of-phase. 


In this application of the amplitude theory, an important experimental and nu- 
merical result is reproduced: the effect of the relative phases of the waves. This is 
an important validation of the approach for this flow which provides encouragement 
for future applications of the theory. Future applications should include detailed 
analysis of the role played by each term in the expansion, and should include three- 
dimensionality as laid out in the appendix. Three- dimensional predictions can be 
compared to three-dimensional numerical simulations as described in §4 below. 


4. Direct numerical simulations 

In the numerical simulations, the Navier-Stokes equations (2.2) are solved using 
a spectral numerical method based on the expansions discussed by Cain, Ferziger 
& Reynolds (1984). Grid resolutions up to 64 3 were used for the results presented 
below. 

In these preliminary studies, we present results for four different simulations. In 
the first two cases (Case I and Case II), the fundamental (a^,0) as well as the 
translational mode (0,7^)are initially nonzero. The fundamental mode is taken to 
be the solution to the Rayleigh equation at a = ap, while the mathematical form 
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FIGURE 5. Evolution of model energies for CASE I, for the modes: o , (a^, 0); 
o , (2o/p,0); a , (0,7f); + (o f , 7 f); “d * (< *f>27f )• Modes (of> 0) and (0,7 f) 
were initially excited with energies 2 X 10 -3 and 1.5 x 10 -4 respectively. 

for the translational mode is chosen to be 

u(*,2/,z,0) =0, 

v(z>y,z y 0) =Ae~ 2y2 cos( 7 fz), 

w(z,y, z,0) =— ^^e~ 2y3 sin(jFz). 

7/ 

The difference in the two cases is the amplitude A . In Case I it is set to 0.025, 
while in Case II it is 0.05. In the latter two cases (Case III and Case IV), the 
fundamental plus the oblique modes (a/p,7j?) and (ajr, —71?) are initially nonzero. 
The normal velocity v for the oblique modes is initially as given above, the other 
velocity components are such that continuity is satisfied and 7 pu — apw = 0. In 
Case III the oblique mode is approximately out-of-phase with the fundamental, i. e., 
it is symmetric about the core of the two-dimensional vortex, while in Case IV the 
oblique wave is in-phase. 

A plot of modal energies versus time for case I is shown in figure 5. In this and 
in the following, the modal energy is defined by 

E(a,y f t) = ^S a>7 J |ui(a,y,7,<)| 2 dy. (3.2) 

Here summation over i is implied, and £»(<*, y, 7, f) is the Fourier transform of 
u i( x yyj z )t) l n the (x,z) plane. The summation sign denotes a sum over the modes 
(«>7)> («,- 7), (-<*,7), and (-a, -7). 




100 


J.J Riley, P. D. Mourad, R . D. Moser and M . M. Rogers 



FIGURE 6. Evolution of model energies for CASE II, for the modes: □ , (aj?,0); 
o , (2aj?,0); * , (0,7f); + («f,7f); and x («f,27f)- Modes (of,0) and (0,7f) 
were initially excited with energies 2 x 10 -3 and 6 x 10 -4 respectively. 

In Figure 5 we see that the energy in the fundamental mode grows, at first expo- 
nentially, then more gradually, and finally levels off at time 12, which corresponds 
approximately to vortex rollup. At the same time the harmonic (2a, 0) grows very 
rapidly, peaking at about t = 9, and then decays as rollup is completed. 

The translative mode, which is only algebraically unstable without the finite 
amplitude, two-dimensional mode, initially grows very rapidly. Further analysis 
of the results showed that this growth is initially mainly in the u-component, and 
is due to the streamwise vortices inducing motion in the y-direction, along the 
gradient of the mean velocity (Squire mode). As the mixing layer rolls up, vortex 
stretching occurs along the layer, resulting in further increase in this translative 
mode. We see that subsequently the oblique mode begins to grow rapidly, and 
ultimately overtakes the translative mode. Analysis of visualizations (see below) 
suggests that this growth is due both to the tilting of the streamwise vortices and 
the distortion of the main spanwise rollers by the the streamwise vortices. 

Figure 6 gives the modal energies versus time for Case II, which has the same 
initial conditions as for Case I except that the translative mode amplitude is in- 
creased by a factor of 2. Comparing Figures 5 and 6, we see that the energies in 
the fundamental mode and its harmonic are essentially the same for the two cases. 
From an energetics point of view the two-dimensional rollup is totally unaffected 
by the three-dimensional motions, even though in this latter case the energy in the 
translative mode is within a factor of two of that in the fundamental mode (at 
about t = 2). Previously Metcalfe et al. (1987) observed that, in the interaction of 
an oblique mode with the fundamental, the fundamental mode played a catalytic 
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FIGURE 7, Iso-level surfaces of enstrophy (a,c) and (b,d) for CASE I (a,b) and 
CASE II (c,d). Solid surfaces represent negative values. Contour levels are 1.5 in 
(a,c), —0.3 and 0.3 in (b) and —0.6 and 0.6 in (d). 


role. Its presence was required for the enhanced growth of the oblique mode, but 
its energy was unchanged by the interaction. The increased energy in the oblique 
mode came from the mean flow. Our results indicate that this catalytic trait holds 
for both the fundamental and its harmonic interacting with the translative mode. 

The total energy in the three-dimensional modes is significantly different for the 
two cases. However, examining the results more closely we find that the total energy 
in a particular three-dimensional mode for Case II is almost exactly four times the 
corresponding value for Case I. Since the initial amplitude for Case II is exactly twice 
that for Case I, this is consistent with the interpretation that the development of 
the three-dimensional dynamics is linear. This result gives some indirect support 
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2 tS/AU 

FIGURE 8. Evolution of model energies for CASE III, for the modes: □ , (a F ,0); 
o , (2a F ,0); a , (0,7f); + (<*f,7f); and x (aF,27F). Modes (aF,0) and (aF,7F) 
were initially excited with energies 2 x 10” 3 and 10 -3 respectively. 


to the model of Corcos & Lin (1984) who, although studying the interaction of the 
fundamental mode with an oblique mode, assumed that the dynamics of the oblique 
mode were linear. 

Perspective plots of iso- value surfaces of both streamwise vorticity (a;*) and the 
norm of the vorticity (u> = y/WiWi) for these two cases are shown in figure 7. Note 
that in the plot of w xy the contour level for Case II was chosen to be twice that 
of the level for Case I. We see the development of the streamwise vorticity along 
the breads, and also these vortices being wrapped into the large-scale vortex core. 
Note that the streamwise vortices appear to be somewhat flat due to the straining 
field of the two-dimensional vortices. Also note the appearance of counter-rotating 
vortices in the core, probably due to the bending of the large-scale vortices by the 
streamwise vortices. It is interesting to note that the streamwise vorticity plots for 
the two cases are almost the same, lending strong support to the interpretation that 
the three-dimensional dynamics are essentially linear. This is surprising in view of 
the large distortions observed in the flow field. The plots of w indicate the distortion 
of the large-scale vortices by the streamwise vortices. 

We next discuss the results for the oblique wave cases: Case III, in which the 
oblique wave disturbance was initially approximately out-of-phase with respect to 
the vortex core, and Case IV, in which it was in-phase. Plots of the modal energy 
versus time for Cases III and IV appear in Figures 8 and 9, respectively. First 
note that in both cases the behavior of the energy in the two-dimensional modes 
is essentially identical to that for Cases I and II. This is further evidence that the 




Three-dimensional forcing of mixing layers 


103 



FIGURE 9. Evolution of model energies for CASE I, for the modes: □ , (a/?,0); 
o , (2a F ,0); a , (0,7 f); + (<*f,7f); and x (a Fy 2j F ). Modes (a F ,0) and (<* f ,7f) 
were initially excited with energies 2 x 10~ 3 and 10” 3 respectively. 

two-dimensional modes are unaffected by the three-dimensional motions, and is con- 
sistent with the catalytic character of the fundamental mode observed by Metcalfe 
et al. (1987). We further note that in both cases the oblique mode (ckfj7f) grows, 
but not as rapidly as the translational modes in the previous cases. Furthermore, 
in these cases the translational mode quickly grows from zero to be the same order 
as the oblique mode. 

Though the energetics of the modes are unaffected by the phases, examination of 
the flow fields reveals that the details of the fields are quite different. Perspective 
plots of (v x and u) for both cases appear in Figure 10. The time selected is t = 10, 
when rollup is almost complete. In Case III there is very little evidence of streamwise 
vorticity in the braid region; it has mainly been rolled into the vortex core. This is 
due to the fact that, initially, the stagnation point in the braids was near the zero 
in the streamwise vorticity. In Case IV streamwise vortices are evident, similar to 
those in Cases I and II. Comparing the plots of a; we see that the flow in the rollup 
regions is very different for the two cases. 

5. Conclusions 

From this preliminary study of the sensitivity of the mixing layer to three- 
dimensional disturbances, we draw a number of tentative conclusions. First it is 
clear that the mixing layer is sensitive to the wavelength of the disturbances (or fre- 
quency in the spatially-developing flow), the relative phases and amplitudes of the 
disturbances, and to the form of the disturbances. Clearly the disturbance needs to 
be well-defined in order to predict its behavior. 
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Figure 10. Iso-level surfaces of enstrophy (a,c) and u> x ( i,d ) for CASE III (a,b) 
and CASE IV (c,d). Solid surfaces represent negative values. Contour levels are 
1.5 in (a,c)> —0.25 and 0.25 in (b) and —0.5 and 0.5 in (d). 


Secondly, it is important to realize that naturally-occurring temporally-evolving 
mixing layers may develop very different three-dimensional disturbances from those 
of spatially-evolving mixing layers. The principal initial three-dimensional distur- 
bances in the spatial layer are fairly narrow-banded streamwise vortices, whose 
development and breakdown appear to be a principal factor in the transition to 
turbulence of the mixing layer. Their initiation appears to be related to features 
of the experimental facility. In the temporal layer it is unlikely that this mode is 
initially excited, and laboratory experiments do not show evidence of these vortices 
in the braid regions. 

Third, at least during vortex rollup, the two-dimensional modes appear to be 
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almost totally unaffected by the three-dimensional disturbances. In addition, the 
three-dimensional disturbances behave in an almost linear fashion, being only af- 
fected by the two-dimensional rollup and not by their own interactions. This is 
true even though their energy levels are large, of the order of the energy in the 
two-dimensional modes. 

Fourth, the translative instability appears to be faster growing than that of the 
oblique modes. This appears to be due to the growth of streamwise velocity fluctua- 
tions (the Squire mode). Fifth, the relative phase of the oblique mode perturbation 
with respect to the fundamental mode is not important in the energetics, but is 
important in the resulting flow pattern. And sixth, the successful preliminary work 
with the amplitude expansions suggests that they may be a useful quantitative tool 
for future research. 

Finally, this preliminary study suggests additional simulations and mathematical 
analysis. For example, it would be useful to more carefully test the amplitude 
expansion technique by comparing its predictions with the results of direct numerical 
simulations for cases, both two- and three-dimensional, that treat the exact same 
initial conditions. If the theory appears to be adequate, then it could be used to 
better understand the interactions occurring in the simulations, and to determine 
what combinations of initial conditions might lead to desired features of the mixing 
layer behavior. With regard to the numerical simulations, there are a large variety 
of initial conditions which should be explored. Further simulations to determine 
the range of possible outcomes would be useful. 
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In this appendix, we develop a theory for wave- wave interactions between oblique 
waves in supercritical shear flows. A more limited set of theoretical equations based 
upon the same procedure is used in the body of the text to study the interactions of 
a fundamental wave and its subharmonic. The more complete equations discussed 
here were developed during the summer workshop, and will be used in future studies. 

The heart of the theory is a set of nonlinear evolution equations which govern the 
temporal history of the amplitudes of interacting waves. The theory is similar to 
Stuart’s energy integral method, and Lorenz’s derivation of his famous amplitude 
equations for convection. As in these other methods, we manipulate the Navier- 
Stokes equations in a formal way, without benefit of asymptotic arguments and 
with severe Fourier, Orr-Sommerfeld and amplitude truncations. Also like them, 
we hope that these truncations are a good approximation of the full series solution, 
and that they lend at least a qualitative insight into the physics of the problem. 
If we had recourse to asymptotically small parameters, a perturbation formalism 
would determine the relative weights of different physical effects. However, since 
the flow in question is supercritical, one cannot create a small parameter out of the 
governing Reynolds number. Neither has the relative energy of the disturbance to 
the mean flow proven useful as a small parameter. The subsequent evolution of the 
waves in the presence of the mean flow is governed by the ability of the waves and 
mean flow to exchange energy with each other, which, even in an initially linear 
regime, is independent of this parameter. With these obvious choices removed for 
supercritical flows there are no readily apparent small parameters that we can use. 
Supercritical flows generally evolve quickly away from any given initial state to a 
relatively remote state - a process which is not usually amenable to asymptotic 
analysis. With the formalism developed here, we can build tools which are simpler 
to use than the full Navier-Stokes equations, as well as lay bare the underlying 
physical processes. It becomes a powerful tool of analysis when used in tandem 
with experimental observations and/or direct numerical simulations. 

The theory starts with the fully nonlinear, nondimensional Navier-Stokes equa- 
tions, in the coordinate system used in the text. Separate the mean flow (a function 
of y and t alone) from the wavelike part (oscillatory in x and z) of the dependent 
variables, with the following notation: 

U(*,y,z,t) = U 0 (y) + U(y,f) + u(z,y,z,t) (Al.l) 

P = P(v,*) (A 1 . 2 ) 

The vector Uo = (i7o,0,0) defines the original flow field which will house the inter- 
acting waves. We have chosen: 


Uo = erf(v^ry) (A1.3) 

where y is nondimensionalized by the initial vorticity thickness 8. This profile is 
selected because it is the self-similar laminar profile for a temporally-evolving mixing 
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layer. The mean pressure field will remain unspecified, as it does not enter into the 
problem. The vector U = (U y 0,W) represents a time-dependent modification to 
the initial mean flow caused by nonlinear feedback from the spatially oscillating 
part of the field. It is calculated as part of the solution, as well as the modified 
mean pressure, given by P. 

Substituting the expansion (Al) into the Navier-Stokes equations one obtains the 
following: 


(u + U)- V{u + U) + u- Vu (A2A) 
(u + U) * V(u) + u • Vv ( A2.2 ) 

(u + U) • V(w +W) + iTVw (A2.3) 
u- Vw (A2A) 

u -Vw (A2.5) 

with: 

A d 2 d 2 d 2 

8x* ^ dy 2 ^ dz 2 ’ 

and the bar over the non-linear terms indicates the average in x and z (but not 
time) of each term. 

A.l Fourier expansions 

The next step is to write the spatially oscillating components of the flow in a 
highly-truncated Fourier series in x and z. This series is truncated at the beginning 
so that the smallest number of Fourier components need be considered. 

The collective nonlinear behavior of the primary waves with wavenumber vectors 
(a, 0 ), (a/ 2 , 0), (a/ 2 , 7 ) an ^ (a, 7) 1S the ma in focus of this study. (Note that, for 
notational simplicity in this Appendix, we are dropping the subscript F, which de- 
notes the fundamental mode in the main body of this paper.) This theory describes 
the evolution of a flow field by studying components of the Fourier transform of the 
velocity fields. The choice of the primary waves is motivated by the observation that 
these are significant waves: with this theory we will only attempt to explain their 
behavior, and can only suggest why these should be the waves of choice. Accord- 
ing to Unear theory, (a,0) (also called the fundamental mode) initially grows most 
quickly, and (a/ 2 , 0 ) is its subharmonic. The subharmonic is significant because 
after the roll-up of the initial vortices (equilibration of the fundamental wave), the 
subharmonic is observed to grow resulting in the pairing of the two-dimensional 
vortices. The pairing process is dependent upon the relative phases of (a,0) and 
(a/2,0), among other things. Observations show that the pairing is followed by, 
and sometimes concurrent with, the development of three-dimensionality in the 


d 1 

v-t + Uou x + v-q-Uo — + P* = — 

v t + U 0 v x - -^-Av + p y = - 
Re 


w t + U 0 w x - —A w +p z = — 
He 


-Lr 

Re 


Ut-—U yy = - 


W ‘-T' W » = - 
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originally two-dimensional flow. Spectrally, this is represented by the inclusion of 
oblique waves (a, 7) and (a/2,7). 

Other waves could also be chosen as primary: for example, an interesting choice is 
the wavenumber vector (0,7), which denotes a system of longitudinal vortices. The 
formalism would proceed in the same fashion if these were included, with different 
evolution equations as the result. However, since this is as yet only a temporal 
theory, and (0,7) seems to be more relevant for spatially developing flow fields, we 
ignore this mode for now. 

There is another class of waves which are included here, because they are forced 
by various combinations of the primary waves and their complex conjugates. For 
short times, these waves are observed to be slaved to the primary ones, that is the 
dynamics of these modes are dominated by the forcing. Their intrinsic dynamics are 
considered insignificant so no extra degrees of freedom will be required to account 
for their effects. For this limited study, we only include the first harmonic of (a, 0), 
given by wavenumber vector (2a, 0) because it is important for the interaction 
between (a,0) and (0,7). We also assume that the modification to the mean and 
the so called Squire modes (see below) remain slaved to the primary waves. 

Which waves to include, and whether they are included as free or forced waves, 
are subtle, contentious, important and open questions. They represent large areas 
of possible improvement as well as additional complexity. 

Finally note that while we have only been discussing modes of the form (p, q), with 
p and q positive, all the modes (±p, ±q) are important. However, the velocity fields 
are real, which implies that two of these modes are not independent due to conjugate 
symmetry. In addition, we will require that the velocity be mirror symmetric in 
the z direction, which will eliminate the independence of another mode. This is a 
symmetry which is preserved by the Navier-Stokes equations so that with an initial 
condition satisfying this symmetry the solution will remain symmetric. Thus with 
the reality of the velocity and the mirror symmetry, only the modes (p, q) with p 
and q positive are independent. 

For each Fourier mode (p, q) the y-dependence of the three velocity components 
must be represented. However, the constraint of continuity eliminates a degree of 
freedom so that only two y functions need be considered. To do this we define two 
new functions of y: 

tp = u? = pw — qu, (.43) 

0 

where u, v and w represent the Fo urier co efficients of the velocity components for 
the given wave number and 6 = yjp 2 + q 2 . This notation is chosen because xp can 
be interpreted as a stream function and is proportional to the y component of the 
vorticity; is also referred to as the Squire mode. Given ip and & for a particular 
Fourier mode, the velocity components can be recovered as follows: 


u = 




qu> 


v — —iSip; w = 


+p u> 


(M) 
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The equations for V 1 and u> for our Fourier mode (p, q) can easily be derived from 
the Navier-Stokes equations, with the result; 


~ + ipUo 

& +ipU ° 


A V> — ipi!> 


(PUo 

dy* 

dU 0 


- -Lav = - 

Re W 6 
1 


d ( i pT x +_igT z ) s2 1 


dy 


u ~ " Te V2 “ = [qF * " Pr * ] <™> 5 


(A5.1) 

(Pt9) 

(>15.2) 


where, T = (U + u) • V(U + u), the subscript on T refers to the vector component, 
A is the Fourier transformed Laplacian (A = S 2 /dy 2 — p 2 — q 2 ) and 8 2 = p 2 + q 2 . 
The notation [*]( p , 9 ) indicates the (p, q) Fourier mode of the quantity in brackets. 
The above equations are valid for all wave number pairs except (0,0). The Fourier 
transformed Navier-Stokes equations for V and W are used directly for this mode. 

We now introduce real amplitudes A(f), B(t ), (7(f), D(t) and phase functions 
0A(t), 0s(t)i 0/?(f)> one f° r each of the primary waves. The ultimate goad is to 

derive ordinary differential equations which govern the evolution of these functions. 
The forced first harmonic of (a,0), the forced modifications to the mean flow due 
to the various waves and the Squire modes do not have independent amplitudes. 
Rather, they are functions of the amplitudes of the forcing waves. This is consistent 
with the equations, and is as described in Herbert (1983). The expansion for u 
is given below. Those for v and w are similar in form. Note that the Fourier 
mode (2a, 0) is especially busy, with two contributions. The functions Ux I (2a l o) an d 
^c,(2«,o) represent the response of the (2a, 0) mode to forcing by the (a,0) mode 
and (a, 7) mode respectively. 


u{x,y,z,t) = 2Real {A(<)u (a>0 )e <(a *-^ ( *» + B{i)u( al2fi) e* a *l 2 - 6 »W 
+ [^ 2 (<)«A,( 2 «,o)e- 2< ^ (<) +C 2 (<)« Ct ( 2 a,o ) e- 2Wc( ‘ ) ] e 2< “'} 

+ 4C(i)Real cos(7z) (^6) 

+ 4£(<)ReaJ {u (a/2t7) e*'(“*/ 2 - <,c ( t ))} cos( 7 z) 

+ A 2 {t)U A + B 2 (t)Us + C 2 {t)U c + D 2 (t)U D , 

where the dependence of u and U on y and t is understood. 

A. 2 Streamfunction equations 

By substituting expansions like (A6) into the Navier-Stokes equations and con- 
verting to and u; as in ( A5) we can write the truncated modal equations governing 
each Fourier mode. In writing these equations we will make use of the following 
operators: 
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£(,„> = i(p‘ + 

+ (A7.2) 

W>(p, q ),ii>(i',p)) = ip [Ap«)-^ M (*,p)A*,p) “ ( M (J ».ff)V»(„,„)] ( 47 . 3 ) 

N(p,q;*,p) = Up + *) (***<«>»<*, ,0 + '.*>) + 

- ^ «)“(*.?) + <W^«(*,p) + *>"(p, 9 )“(»,p)) ( A7 - 4 ) 

Below are examples of the resulting equations governing the stream- functions 
Equations for the Squire modes w are similar. Note that the non-linear terms are 
expressed (using the above operators) in terms of the stream function whenever pos- 
sible, but that the terms involving the N operator require the velocity components, 
which are found using (A 4 ). Also note that these are partial differential equations 
in y and t . First is the equation for the stream function of the fundamental wave. 

( 3 \ 2 

ft + o(<) j M( a(0 )V>(cr,0) + ■£(a,0)V’(a,0) = -^e - ‘ (2 * s ~° A ^(^ a /2,0), V*( a/2,0 )) 

+ [JV (a/ 2 , 7; a/ 2 , —7) + N(a/ 2 ,- 7 ;a/ 2 , 7 )] 

+ A 2 [{ft(^(a,o),^,(2cr,0)) + ^A,(2a,0), V»(*a,0)) } + O(V>(a,0) , ^X,(0,0))] 

+ C 2 [|ft(V , (* a ,0)>V , C,(2a,0)) + ^(V , C,(2a,0)> V’fa.O))} + V’C,(0,0))] 

+ [S 2 n(V>(a,0),V’B,(0,0)) + £> 2 ft(V>(or,0),V’D,(0,0))] * 

. (A8,1) 

Where ^(o.o) = / U(y')dy' and the subscript letter on ^(0,0) and ^(20,0) indicate 
which primary mode is providing the forcing. The five different terms on the right- 
hand side of (A8.1) represent the following effects on the fundamental (in order): 1) 
the effects of the subharmonic 2 ) the influence of the oblique subharmonic, 3 ) the 
combined influence of those parts of the forced first harmonic and the modification 
of the mean flow which are caused by the fundamental wave, 4 ) the combined 
influence of those parts of the first harmonic and the modification of the mean flow 
which are forced by (a, 7), and 5) the influence of those modifications to the mean 
flow which are induced by the other waves. Next are the equations for the two parts 
of the forced first harmonic, 

+ ^(2a,0)^A,(2a,0) + £(2a,0)^A,(2a l 0) = ^(^(a.O) > ^(a.O)) (A8.2) 

(Jj: + 2c(<) ) M (2°.0)V’C',(2a,0) + L(2c,0)‘<l>C,(2a,0) = W(a,7;a,-7) + N(a, —7; a, 7) 

(48.3) 
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The first equation (A8.2) describes the part forced by the fundamental (a,0), while 
the second equation (A8.3) describes the part forced by the oblique fundamental 
(a, 7). Finally, we present the equation for the modification of the mean flow by 
the fundamental mode. 


(I + w Ua + l ™ Ua = b *(...>)] 

(.48.4) 

For equations (A8.2,A8.3 and A8.4) we are only interested in particular solutions 
since these are modes which are slaved to the primary modes. Equations for modes 
not presented above and equations for u> have similar form. 

Functions a, 6, c and d which appear in the above equations are defined below. 
They are the left-hand sides of what will become coupled Landau equations, which 
will give temporal development of the amplitudes and phases of the disturbances. 

< A9 - 2 > 

< 49 ' 3 > 

< MA) 


A.2 The amplitude expansion 

The next step in the development of the theory is to develop expansions to de- 
scribe the time and y dependence of the V>’s and the time dependence of the functions 
a, b , c and d. In choosing these expansions we arrange the terms so that the depen- 
dence on the amplitude and phase functions A, B, C , D, 6 a, $b, 6c and 6d, of the 
left and right hand sides of equations like (A8.1) will match. Thus the expansion 
for V>(a,o) is 


B 2 


= ^O.(a.O) + * a) + 4 > 2,(a,0)-J 

+ [A 2 ^3 t ( Q , 0 ) + B 2 4>4,(a,0) + C 2< l>S,(a,0) + D* <f>6,(a,0)] 

+ [<t>Ua,o)e i(0A - 9B ~ ec+6D) +<t>s,( a ,o)e i{eA+<>B - 0c+eD) 

+<f>9,(a,O)e i(eA - 3dB - 9c+0D) ]+..., 


£l e -H20 D -e A ) 


(All) 


where it is understood that the <f>j are functions of y. Note that the amplitude and 
phase dependence of the first several terms is the same as the coefficients appearing 
on the right hand side of (A8.1). In addition there are the terms with coefficient 
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BCD/ A y these terms arise because of the time derivative operator on the left side 
of (A8.1), that is the term 


(Jt + o(<) ) 

If we had used the ‘shape’ assumption of Stuart or Lorenz, (All.l) would consist of 
only the first term in this expansion, which is the solution to the Orr-Sommerfeld 
equation (see below). The current expansions retain more of the relevant physical 
processes. Similarly we obtain expansions for the time functions which have ampli- 
tude and phase dependencies dictated by the form of the equations like (A8.1). For 
example compare (All) above with (A12.1) below. In the following equations the 
coefficients a», 6*, c», and di (Landau constants) are complex constants. 


a(t) = do + aj-j-e + <12 — j-e + [<13 A 2 + a^B 2 + a$C 2 + ae-D 2 ] 


£< 17 ~**c+ Od) _j_ -0c+Qd) 




(A12.1) 


b(t) = bo + biAe-* 9 *- 20 *) + b2^-e- i(9c - 0D - 9B) + [ b 3 A 2 + b 4 B 2 + b 5 C 2 + b*D 2 ] 

B 

+ b 7 A 2 e- 2 * 0A - 29o) + b t D 2 e- 2i{ ° D - 0B) 

+ A 2 e ~ , ( 0A ~ 20B ) [69 cos(0a — 2 Ob) + 6, 0 sin^x — 26b)] 

+ A ^ D ^b n e- i(0A - 30 B+Oc-e D ) + bl2e i(0A+e B -0c-9 D ) 

+ b 13 e~ i(9A - 0B - e ° +9D) + e - ii9A ~ 29B) {b 14 cos(0b - 6 C + 6 D ) 

+6,5 sin(<? B - Oc + Od)}] + • • • (^112-2) 

TIT) 

c(f) = c 0 + ciA 2 e~ 2 ^ 0A ~ 0B ) + c 2 —p^-e~ l( ' 0B+9l> ~ 9c ^ + [ c 3 A 2 + c 4 B 2 + c 5 C 2 + c 6 D 2 ] 

c 


_ l — ^ L _ ^z-Wa-Qb-Oc+Od) c%e -i(o A ^e B -e c ~e D ) j + ^ 


(A12.3) 


d{t) = d 0 + d^Ae-^- 20 ^ + d 2 ?f/e-* 6c - 0B - 0D) + [ d 3 A 2 + d 4 B 2 + d 5 C 2 + d 6 D 2 ] 


+ d 7 A 2 e- 2i(0A ~ 20D) + d 3 B 2 e~ 2i(0B - 0D) + d { 

ABC 


ACD 

-B~' 


} —i{0 A —3 0 B + 0 c -~&d ) 


_j_ ^dioe~ i ^ A ~* B + $ c- 3$d) ^^^(Oji-Ob-Oc+Od) 

+d 12 e - i(0A+0B - 0c - 0D) ] + ... 


(A12.4) 


Combining equations (A9) and (A12), ordinary differential equations for the am- 
plitude and phase functions are obtained. It is these evolution equations that must 
be solved to obtain the time dependence of the desired solution. All that remains 
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is to compute the evolution coefficients a* etc. The methodology for computing 
them is discussed below; however, they have not yet been computed for this three- 
dimensional case. Still, we can make some qualitative inferences from the form of 
these equations. Consider (12), in the form of Aa(t ), Bb(t), etc., so that one can 
consider dA/dt ... If >1(0) ^ 0 while B( 0), C(0), D( 0) = 0, only the fundamental 
mode will evolve, since the other modes will remain zero. Their initially zero growth 
rate ( dB/dt(t = 0) = 0, etc.) will remain zero ( dB/dt = 0, etc. Vt). Eventually, 
the fundamental will equilibrate due to the feedback from the modification to the 
mean flow and forced first harmonic, embodied in a3i4 3 . This agrees with obser- 
vations and the numerical simulations of Metcalfe, et al . (1988) among others. If 
only the fundamental and subharmonic have nonzero initial amplitudes, they will 
interact, causing pairing; however, the three-dimensional modes will remain zero, 
again because their forcing requires that at least one of them be nonzero in the 
initial conditions to get a nonzero time rate of change. Consider the case of nonzero 
initial amplitudes, save for the oblique subharmonic (a/2,7). The fundamental 
and subharmonic will interact as before, while the more slowly growing, or even 
decaying oblique fundamental (<*,7) is initially of little consequence. The oblique 
subharmonic will remain negligible, because it requires the combined effort of the 
two-dimensional subharmonic and oblique fundamental to stimulate it into growth, 
represented mathematically by the d 2 term in (A12.4). When the fundamental gets 
large enough, its forced first harmonic will stimulate the oblique fundamental into 
accelerated growth (represented mathematically by the C\ term in (A 12.3)). This 
should occur at least simultaneously with, if not after pairing, because the a\ and b\ 
expressions which govern pairing are larger than the C\ term, assuming comparable 
sizes of these coefficients. Finally, after the oblique fundamental mode reaches a 
size of some consequence, the oblique subharmonic will appear, driven by the two- 
dimensional subharmonic and the oblique fundamental via at least the d 2 term. It 
is too difficult to judge what happens beyond this. This scenario, suggested by the 
equations, is similar to numerical simulations and observations. Actual integration 
of the evolution equations is required to verify these conjectures. 

One can test the individual effects and hence the importance of each of these terms 
through a series of integrations of the equations with and without the expression of 
interest. This will allow the analysis of the roles played by the various interactions. 
Comparisons with Corcos and Lin (1984) would be interesting, since they considered 
a partially linearized form of the Navier-Stokes equations, which contain, some, but 
not all of the interactions built into our equations. Future work will consider such 
a program of experimental mathematics. 

Evaluating the expansion coefficients 

Substituting expansions like (All) and (A12) in (A8) and equating like terms 
in the amplitude and phase functions results in a series of linear equations for the 
expansion functions <j>i . The equations for 4>o are Orr-Sommerfeld equations. For 
example, the (a,0) equation is 


aoAf( Qt o)<£o,(a,0) + £(a,0)<£o,(a t 0) = 


(A13.1) 
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where a<> is the Orr-Sommerfeld eigenvalue and <f> o is the eigenfunction. We will 
take as the solution the first (most unstable) eigenfunction and eigenvalue. The 
rest of the equations are linear, inhomogeneous and of the form 

(Vj + <*o)M(a,0)^i,(a,0) + L(a,0)4>j,(a,0) = Yj - OjAf( a)0 )^ 0| ( a>0 ), (i413.2) 

where the rjj are nonzero complex numbers, and the Yj are functions of y. Both rjj 
and Yj can be determined from the solutions already obtained for a*, b{ etc. and 
<j>i for i < j. As an example, the remainder of the analysis will be done for the 
fundamental mode (a,0); the analysis for the other modes is similar. 

Since the functions are only defined to within a multiplicative constant, a nor- 
malization condition is necessary to obtain a solution to these equations. Also, since 
the linear solution ^o,(a,o)(y) has an arbitrary magnitude, so will the corrections 
( <j>j ). It is convenient to choose a linear norm, and apply it in a way which allows 
the determination of the Landau coefficients. Herbert (1983) chooses a linear nor- 
malization, (designated by {•}), such that the following normalization condition is 
satisfied: 

{V>(a,o )(2/>f)} = 1 for all time. 

We use {V>(a,o)(sM)} = ^(«,o)(y»0^y> although any linear norm which is not 

zero for ^ and <f > o would do. Note that this condition is not simply {<£o,(a,o)} = 1. 
This unusual choice of normalization places conditions on all of the <£ nt ( 0f o)> namely: 

{4 W)(SM)} = {E ^n(W«,(c,o)(y)} = E F n (t) {4> n , (a ,o)(y)} = i (-414) 

\ n=0 J n=0 

Since the F n (t) are functions of time (dependent on the amplitude and phase func- 
tions, see equation All) and Fq = 1, this condition is satisfied by requiring 

«*.,,}={£ 

Setting the right side of (A13.2) to zero produces the Orr-Sommerfeld equation 
with a parameter different from the first eigenvalue ( r)j / 0); therefore, equation 
(A13.2) does not have homogeneous solutions (unless rjj + ao happens to also be an 
eigenvalue). The solution to (A13.2) is therefore just the sum of the responses to 
the individual inhomogeneous terms: 


= Tj - ~~4 > 0,(a,0)i (>116) 

VJ 

where Tj is the solution to Equation (A13) corresponding to Yj , and aj/r)j<f> o,( a ,o) 
is the solution corresponding to the term on the right side of (A13.2). Applying 
the normalization condition to (A16) and using (A15) the value for the Landau 
constant aj is obtained, 

a j = {Tj} 7 }j* 


0*18) 
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Thus we have found both the Landau constants and the expansion functions 
<f > . In principle this procedure could be carried out to infinity, but in practice a 
severe truncation would again be applied. With the Landau constants, the ordinary 
differential equations for the amplitudes and phases can be integrated, and used to 
construct the velocity field in space and time. 



